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Abstract 

We review the definition of the Casimir energy steming naturally from the con- 
cept of functional determinant through the zeta function prescription. This is done 
by considering the theory at finite temperature and by defining then the Casimir 
energy as its energy in the limit T — > 0. The ambiguity in the coefficient (7^/2 is 
understood to be a result of the necessary renormalization of the free energy of the 
system. Then, as an exact, explicit example never calculated before, the Casimir 
energy for a massive scalar field living in a general (H-2)-dimensional toroidal space- 
time (i.e., a general surface of genus one) with flat spatial geometry — parametrized 
by the corresponding Tcichmiiller parameters — and its precise dependence on these 
parameters and on the mass of the field is obtained under the form of an analytic 
function. 
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1 Introduction 



The Casimir effect ||T]] is a beautiful and simple manifestation of the influence that bound- 
aries or non-trivial spacetime topologies have on quantum field theories (see for example 
0, ^, ^, H). The modern approach for the calculation of the Casimir energy is the zeta 
function regularization scheme put forward in The idea is the following. In order to 
have a well-defined notion of energy, let us work in a (i-dimensional ultrastatic spacetime 
0. Thus the metric in some coordinate system is of the form g = —{dx^)'^ + gs with 
the spatial part gs of the metric (see Eq. ( p.l|) ). The differential operator D describing 
the field equation may be decomposed as D = —Oq + Ds- Introducing E"^ = as the 
eigenfrequencies of Z^^, the zero-point energy is formally given by 

^ n 

It may be regularized by defining 

Ereg{e) = ^fi'Xsi-l/2 + e), (1.2) 

with Cs{s) being the zeta function associated with the {d — l)-dimensional operator Ds- 
The scale /i with dimension (length)"^ has to be introduced in order to keep the zeta 
function dimensionless for all s. 

General zeta function theory p[ tells us, that Ereg{^) is a meromorphic function with 
a pole at e = 0, its residue being — (l/2)Cd/2(-Ds')/(47r)'^/^. The coefficient Cd/2{Ds) is 
the Seeley-De Witt coefficient appearing in the asymptotic expansion for small t of the 
heat-kernel associated with Ds, 



m = E 



e 



-\nt 



^ E CiiDsY. (1.3) 



KA-KtJ 



/=0,1/2,1, 



The pole appearing in Eq. ( |1.2| ) has to be absorbed into the bare action which thus must 
contain a term proportional to Cd/2{Ds)- It is clear then, that the Casimir energy has an 
ambiguity proportional to Cd/2{Ds) (we will come back to this point later). 
Adopting the minimal subtraction scheme, one defines 



Ecasim^r = J llm J/i^^ [C5(-l/2 + e) + C5(-l/2 " e)] 



;i.4) 



where the symbol PP stands for taking the principal part. Based on this definition, 
during the last years the Casimir energy has been calculated for a variety of examples 

i, ^M- 
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The definition ( |1.4| ) of tlie Casimir energy is completely equivalent to tlie one steming 
naturally from the definition of functional determinant by the zeta function prescription 
(for very recent considerations on this issue see [111]). This may be done by considering 
the theory at finite temperature and by defining the Casimir energy as its energy in the 
limit T ^ 0, this idea going back to Gibbons [|12| who considered the single quantum 
mechanical oscillator. In the more general context of quantum field theory under some 
external conditions like boundaries or gravitational fields this definition has been employed 
for example in [0, |1^, |T5|, |16[)- The ambiguity in the coefficient (7^/2 may be understood 



to be a result of the necessary renormalization of the free energy of the system. 

Having summarized the main arguments in favor of the definition of the Casimir energy 
as given by Eq. (|1.4|) , in section 3 we will present the calculation of the Casimir energy 
for a massive scalar field in a general (1 + 2) -dimensional toroidal spacetime with flat 
spatial geometry. The general fiat geometry will be parametrized by the corresponding 
two Teichmiiller parameters and the complete dependence of the Casimir energy on these 
parameters and on the mass of the field will be obtained under the form of an analytic 
function, by using the extended Chowla-Selberg zeta function formula derived by one of 
us in Refs. [p]7| , p^ . For the massless case we obtain complete agreement with previous 
results by Dowker |19]. 



2 Zeta function definition of the Casimir energy 



Let us first briefly summarize the motivation for the deflnition ( |1.4D of the Casimir energy. 
For deflniteness, let us consider the quantum fleld theory of a free scalar field in curved 
spacetime, the Dirac field may be treated analogously. As mentioned before, in order to 
have a well defined notion of energy we shall restrict our considerations to a cZ-dimensional 
ultrastatic spacetime Ai, possibly with a boundary, and with the metric 

ds"^ = -dxl + gab{x)dx''dx'', (2.1) 

where x = (xi, ...,Xd-i)- The action of the field theory we consider is [0], pO| , 



S = -ll d'xlgl^^x) (□ - ^i? - m') 0(x), (2.2) 

M 

with □ being the Laplace-Beltrami operator of the ultrastatic spacetime. Variation of 
equation ( |2.2| ) subject to the constraints 

6(j){x) = 0, 
n^'V^'6(j){x') = 0, 



where the prime refers to quantities defined on the boundary dAi, yields the equation of 
motion 

(n-^R-m^)^{x) = 0, (2.3) 
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which is the generahzed Klein-Gordon equation. The following discussion will be quite 
general, so the boundary condition need not to be specified at this point. A unique 
boundary value problem is posed, for example, by assuming Dirichlet- or Robin-boundary 
conditions on the field. 

In an ultrastatic spacetime the quantum field theory at finite temperature may be 
developed in complete analogy with the Minkowski-space theory and that is why we will 
skip the details of the calculation. In the Euclidean formulation of the finite temperature 
theory, the partition sum Z may be written under the form of a functional integral of the 
kind 

Z\i5\ = |[rfy.]exp|-i(^,D^)}, (2.4) 
where we have used the scalar product 

r 

{f,h) = JdTjdJ:\g\'/'rh (2.5) 
s 

for the two vectors / and h, and being S the spatial section of the manifold Ai. Here 
the integration extends over all fields periodic in the imaginary time r with periodicity 
P = 1/T and fulfilling the boundary conditions at the spatial boundary. The operator D 
is given by 

D = -—- A + ^R + m\ (2.6) 

Then, one formally obtains 

Z[l3,ij] = (2.7) 

the scale A |^ being necessary in order to keep everything dimensionless. The functional 
determinant of the operator X^^D needs, of course, regularization. We will use the zeta- 



function regularization scheme introduced by Dowker, Critchley |^ and Hawking 
In this scheme, the free energy is defined as 

m = [O(0,/3)lnA2 + C(0,/5)] . (2.8) 

The function (d{s,P) is the zeta-function associated with the operator D. Using the 
ansatz 

1 

= - exp 

the eigenvalues are seen to have the form 




— +El leTL, 



P J 

with the energy eigenvalues being defined through 

(-A + ^R + m^)M^) = ElM^)- (2.9) 
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Making use of a Mellin-transformation and a theta-function identity the free energy 
may be written in the form 

= lppCs(-l)+-^CJ\nX'-l + 2\n2] (2.10) 



2 V 2 



2(47r) 



e 



The energy of the system is then given by 



E=^PF[P]. (2.11) 



Finally, defining the Casimir-energy as the limit of (|2.11|) for T ^ 0, we obtain 



Ec^lSM] = lin, £ = ipPC, (-i) + 1„A=, (2.12) 

with A = We thus arrive to the definition of Refs. [0, |^. In the last reference a 
detailed discussion of the meaning of this definition and of the problem of renormalization 
has been carried out. From the above derivation of the definition of the Casimir energy 
it is completely clear that the ambiguity of the Casimir energy is simply a result of the 
(in general) necessary renormalization of the free energy. In the case when the coefficient 
Cd/2 vanishes, the definition ( p. 121) gives a well defined, finite value. 



3 Casimir energy in a (l+2)-dimensional toroidal space- 
time 

Let us consider, as an example, a (H-2)-dimensional spacetime with the topology R x 
2M. We will concentrate on the case when the geometry of the space S ~ is locally 



flat. One can construct such geometry, as is usual, by taking E = [0, 1] x [0, 1]/ ~, 
where the equivalence relation is defined by (^i,0) ~ (^i,l) and (0,^2) ~ (1;^2)- A fiat 
2-geometry is endowed on E by giving it a metric 

ds^ = KbdCde, (3.1) 

where 

= |;|,). (3.2) 

The (ti, T2) are the Teichmiiller parameters, independent of the spatial coordinates (^1, ^2); 



and r = Ti + 2x2, T2 > [25, 26 



The Laplace-Beltrami operator of this metric is given by 

1 



A = {\T\'dt -2ndid2 + d^), (3.3) 

r2 
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being its eigenvalues 

V,n2 = (kP'^i - 2x1^1^2 + nl). (3.4) 

In the massive case, m ^ the spectrum runs over ni,n2 € TL. In the massless case the 
zero-mode of A, rii = n2 = 0, has to be excluded. 

An exact analysis of the Casimir energy for this spacetime is possible since it reduces 
to a case of the Chowla-Selberg zeta function (when m = 0) or to one of the extended 
formula that has been obtained recently (case m 7^ 0). 

In fact, in the case when m 7^ 0, the relevant formula is a particular application of the 
following. Let us consider the double series 

E{s]a,b,c;q)= ^ {am^ + bmn + cn^ + q)~^ , (3.5) 

with g 7^ (in general), the parenthesis in ( p.5|) is the inhomogeneous quadratic form 

Q{x,y)+q, Q{x,y) = ax"^ + bxy + cy'^, (3.6) 

restricted to the integers. In the general theory that deals with the homogeneous case, 
one assumes that a, c > and that the discriminant 

A = 4ac - 6^ > (3.7) 

(see [0). Here we will impose the additional condition that q be such that Q{m^n) + 
g 7^ 0, for all m, ra G Z. In the usual applications of the theory, those conditions are 
indeed satisfied. For the analytical continuation of ( p.5|) , the following expression has 
been obtained in Refs. O, IT8[ 



E{s- a, 6, c; q) = 2Ceh{s, q/a) a'^ + ^'^^-1/2 ^(s - 1/2)Ce//(s - 1/2, 4ag/A) (3. 



+— E ^ ' cos{n7rb/a) E ^ ^ + — \/ A + 



J2d\n denoting the sum over the divisors of n and where the function Ceh{s,p) (the one 
dimensional Epstein-Hurwitz or inhomogeneous Epstein series) is given by 

Ceh{s;p) = E(^'+P)"' = ^E'(^'+P)"' (3-9) 

Eq. ( |3.(j| ) provides the analytical continuation of the inhomogenous Epstein series, in the 
variable s, as a meromorphic function in the complex plane. Its pole structure is explicitly 
given in terms of the well-known pole structure of Ceh{s,p)- 
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Eq. ( p.8| ) has been found by one of us and called the extended Chowla-Selherg formula, 
since it contains the Chowla-Selberg formula as the particular case (? = 0, i.e., 



E{s- a, b, c- 0) = 2C(2.) a"^ + ^^^^^ - l/2)C(2s - 1) + 

X 5]^n^-^/Vi_2.(n) cos(n7r6/a) i^s-i/2 (^^^j • (3. 10) 



where 



Formula ( |3.8| ) has been obtained for the first time in , and a misprint has been corrected 
[1^. The good convergence properties of expression ( p.lO|) , that were so much prised 



a,{n)=Y,d\ (3.11) 

d\n 



m 



by Chowla and Selberg, are shared by its non-trivial extension ( |3.8| ). This renders the use 
of the formula quite simple. In fact, the two first terms are just nice — under the form 
( ^.9| ) — while the last one (impressive in appearence) is even more quickly convergent 
than in the case of Eq. ( p.lQ ), and thus absolutely harmless in practice. Only a few first 
terms of the three series of Bessel functions in (|3.8| ), ( ^^ need to be calculated, even if 
one demands good accuracy. We should also notice that the only pole of ( |3.10| ) at s = 1 
appears through C{2s — 1) in the second term, while for s = 1/2, the apparent singularities 
of the first and second terms cancel each other and no pole is formed. Analogously, the 
pole at s = 1/2 in (|3.8| ) comes only from the first term. Eq. ( p.8|) also has these good 
properties, for any non-negative value of q. In fact, for large q the convergence properties 
of the series of Bessel functions are clearly enhanced, while for q small we get back to the 
case of Chowla and Selberg. Notice, however, that this is not obtained through the high-g 
expansion (e.g., just putting g = in (|3.^ )), but using a low-g, binomial expansion of the 
kind 



[a{n + of + q]-' = a-^ ± ^ Xf^"^^'^ {-Y Ch{2s + 2m, c), (3.12) 



which is convergent for q/a< 1. For g — it reduces to a~*C//(2s, c). Actually, formula 
( p.8|) is still valid in a domain of negative g's, namely for q > — min(a, c, a — b + c) . 

Turning now to the particular application of the formula in our specific situation, we 
see that the zeta function corresponding to the Laplace-Beltrami operator in the massive 
(resp. massless) case is simply given by: 



— s 



and 



CA+m<s) = m-'^ + l^] E U\r\^ -2n,l;^) (3.13) 



Ca{s) = (^) E [s- |r|2, -2ri, 1; o) , (3.14) 
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respectively. The values at s = —1/2 are finite and define the corresponding Casimir 
energy. After performing the necessary calculations, and according to the prescription 
that has been derived in the first part of this paper, we obtain as result the quite simple 
expressions 

(3.15) 



n=l V TT . 



-8x2 "^TL ^ cos(27mx^) E d'^ 



1 + — — Ki 2-nnx 



m? 



4^ r (2vrx2rf)'^ f^^^V (2-X2rf) 
and 

CO 

Ca(-1/2) = -— + 47rC'(-2)x^ - 8x2 E cos(27rnx^)a2(n) Ki (2™^) , (3.16) 
in terms of the variables 



x, = ^^, -2 = ^^. (3.17) 

The extrema of the corresponding Casimir energy for the case = 100, in terms 
of the original Teichmiiller coefficients ri and T2, are to be read from Figs. 1-3. In the 
three-dimensional plot over the plane ri, T2 (Fig. 1), the maximal Casimir energy is seen 
to be localized on the section T2 = 1. In order to show this fact more clearly, in Fig. 
2 we have represented the section ri = 0, but the situation is common to any section 
Ti = const. On the section r2 = 1 a periodic structure appears asociated with the value 
of Ti along this section (Fig. 3). This behavior is easy to recognize from the form of 
the function (^{—1/2), (|3.16|) , and is common to any section T2 = const. All the figures 
depicted here have been obtained taking the first 20 terms from the sum over n in ([3.161). 
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Figure captions 
Fig. 1 

Three-dimensional plot over the plane ri, T2, of the Casimir energy for the case m = 10. 
One clearly observes that the extremal Casimir energy is localized on the section T2 = 1. 
All the figures have been obtained taking the first 20 terms from the sum over n in ( p.l6| ). 

Fig. 2 

Plot of the section ti = 0, where the fact that the extremal Casimir energy is obtained 
for r2 = 1 is seen undoubtedly. This behavior is common to any section ri = const. 

Fig. 3 

Plot of the section T2 = I. Here a periodic structure appears asociated with the value 
of Ti along this section. This behavior is common to any section T2 = const. 
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Casimir Energy, m=10 
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